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Abstract. We review and clarify the assumptions of our basic model for neutrino produc- 
tion in the cores of quasars, as well as those modifications to the model made subsequently 
by other workers. We also present a revised estimate of the neutrino background flux and 
spectrum obtained using more recent empirical studies of quasars and their evolution. We 
compare our results with other theoretical calculations and experimental upper limits on 
the AGN neutrino background flux. We also estimate possible neutrino fluxes from the 
jets of blazars detected recently by the EGRET experiment on the Compton Gamma Ray 
Observatory. We discuss the theoretical implications of these estimates. 

1. Introduction 

Quasars are the most powerful emitters of radiation in the known universe. 
These remarkable objects are presumably fueled by the gravitational energy 
of matter falling onto a supermassive black hole at center of the quasar core, 
though the mechanism responsible for the efficient conversion of gravitation- 
al to observed luminous energy is not yet known. If this conversion occurs 
through the acceleration of high energy cosmic rays, their interactions can 
lead to the copious production of hadronic and subsequent leptonic high 
energy byproducts such as neutrinos. It is in this context that quasars have 
long been considered as potential sources of high energy cosmic ray neutri- 
nos (Berezinsky 1977; Eichler 1979; Silberberg and Shapiro 1979; Berezinsky 
and Ginzburg 1981). 

In related work, Biermann and Strittmatter (1987) suggested that pro- 
tons could be accelerated to ultrahigh energies in the hot spots found in the 
jets of radio galaxies (see also Rachen, Stanev and Biermann 1993). Subse- 
quent X-ray production at these sites was considered by Mannheim, Kruells 
and Biermann (1991). High energy 7-ray and neutrino production in jets was 
explored by Mannheim, Stanev and Biermann (1992), following the impor- 
tant discovery of strong high energy 7-ray emission from the blazar 3C279 
(Hartman, et al. 1992). New results on the detection of large fluxes of high 
energy 7-rays from ~ 40 such sources, obtained using the EGRET detector 
aboard the Compton Gamma Ray Observatory (Fichtel, et al. 1994), have 
made the study of these high energy sources even more interesting. 

An important step in modeling quasar neutrino production process was 
made when Begelman, Rudak and Sikora (1990) and Stecker, et al. (1991) 
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pointed out the importance of photomeson production off thermal UV accre- 
tion disk photons as the most significant production channel. 

In our work (Stecker, et al. 1991), we pointed out that following photopion 
production, pair cascading off the intense quasar photon fields would both 
destroy accompanying high energy 7-rays and produce the quasar's X-ray 
emission (e.^.. Done and Fabian 1989). Thus, by normalizing to the X-ray 
emission, we argued that active galaxies could be by far the strongest sources 
of ultrahigh energy neutrinos in the universe, producing fiuxes detectable 
on the Earth with present technology. Measurement of these neutrino fiux- 
es and spectra could thus provide key information as to the fundamental 
astrophysical nature of quasars. 

We further pointed out that by integrating over all quasars in the Uni- 
verse, a large and measurable isotropic neutrino background would be pro- 
duced. Because of the degradation of 7-rays to X-ray energies, the energy 
in quasar background neutrinos could be comparable to that of the cosmic 
X-ray background rather than the much smaller energy of the cosmic 7- 
ray background. We also noted that radio quiet quasars are a factor of 10 
more abundant than their radio loud counterparts (Sramek and Weedman 
1978; Kellerman, et al. 1989). They thus take on corresponding significance 
when considering their contribution to the high energy cosmic neutrino back- 
ground. 

Szabo and Protheroe (1994) have followed up on this basic model with 
independent calculations. They have examined a series of models and assumed 
that nucleons are trapped indefinitely in the quasar core. In such models, the 
nucleons eventually lose energy and interact with the X-rays in the core as 
well as the UV photons, even though the X-ray field of the core is optically 
thin to nucleons. The result is a much higher fiux of lower energy neutri- 
nos than we calculated. The difference between the predicted event rates 
in neutrino detectors between the two calculations is not as pronounced. 
This is because the detectability of high energy neutrinos (cx cTj/JvX range of 
produced muon) is cx £'j,ln(£'j,) in this energy range. 

An extensive discussion of quasar neutrino production took place at an 
historic workshop in Hawaii, with almost all of the above-mentioned work- 
ers taking part. The reader is referred to the proceedings of this workshop 
(Stenger, et al. 1992) for more background into this subject. In that work- 
shop, Stecker, et al. , Szabo and Protheroe, Biermann, and Sikora and Begel- 
man all gave estimates of the neutrino background from AGN. 

In this paper, we will update our model using new information on quasar 
evolution and will also consider neutrino production in blazar jets. We will 
then compare our results with other work. We will also review the present 
experimental upper limits on cosmic high energy neutrino fiuxes and their 
implications. 
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2. The AGN Phenomenon 

The standard "AGN paradigm" assumes that infalling matter forms an 
accretion disk around a black hole. It has been suggested by Kazanas and 
Ellison (1986) that an accretion disk shock will form at a distance of several 
Schwartzschild radii [R^ = 2GM/c^) from the black hole. First-order Fermi 
acceleration of charged particles at this shock then converts a significant frac- 
tion of the gravitational energy into highly relativistic particles which have 
the characteristic hard power law spectrum typical of shock acceleration. 
Other possible energy generation mechanisms involve electromagnetic accel- 
eration by black hole accretion disk dynamo models. In any model where 
high energy cosmic rays are produced at the core of the quasar, such cosmic 
rays (even neutrons) will be trapped by the intense radiation fields. Thus, 
a necessary consequence of any acceleration mechanism in which protons 
reach high energies is the production of high energy neutrinos through the 
collisions of these particles with the intense photon fields in quasars. 

Recent observational studies bearing on quasar evolution and theoretical 
models of quasar formation and evolution can be used to guide one in choos- 
ing the basic physical parameters needed to calculate neutrino fiuxes from 
individual quasars and also to calculate the neutrino background from all 
"unresolved" quasars. There is empirical justification for assuming that all 
quasars run at some constant high efficiency. Quasar evolution appears to 
be characterised a phase of emission near the maximum allowed rate. This is 
quantified in terms of the the Eddington limit (Ledd = '^irGMm.pc/ aj where 
M is the mass of the black hole), the maximum steady state luminosity that 
can be produced before radiation pressure disrupts the accretion fiow. Stud- 
ies have shown that quasars and Seyfert galaxies are emitting at a substantial 
fraction of the Eddington luminosity at all observed redshifts. For quasars, 
this fraction, which we will call the "Eddington fraction", appears to be 
usually above 10~^, and for the less luminous Seyfert galaxies, above 10~^ 
with a substantial number above IQ-^ (Pad ovanni and Rafanelli 1988). 

These results are in accord with theoretical beliefs that as long as there is 
an adequate supply of gas available, the luminosity produced by the "central 
engine" of the quasar should approach Ledd- This also implies that the bright 
phase of a quasar should not he substantially longer than the Eddington time, 
igjj = Mc^/Ledd = 4.4 X 10® yr. Thus, it appears that continuous models 
of quasar evolution are ruled out (Cavaliere and Padovani 1988). 

Haehnelt and Rees (1993) have pointed out that the similarity between 
the timescale t^dd and the evolutionary timescale for the hierarchial buildup 
of galaxy size masses at high redshifts allows for a natural explanation of the 
strong evolution of the quasar population. In their work, they point out that 
if the characteristicly large Eddington fraction of quasars is insensitive to 
redshift, the only way to understand the rapid increase of quasar luminosity 
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with redshift is to allow that the black holes which are active at low redshifts 
are different sources with typically smaller masses than those which were 
emitting at higher redshifts. 

We thus arrive at a picture where all black holes in their active phase 
had an adequate supply of gas to accrete near their Eddington limit so that 
the Eddington fraction is both large and roughly constant with redshift, as 
we assumed in our original paper (Stecker, et al. 1991). In this scenario, 
typically the larger mass black holes turned on and used up their supply of 
gas early, to be replaced in later generations by smaller mass black holes, 
which, in turn, used up their gas. Such a picture would naturally explain why 
the quasar X-ray spectral index is no different at 2: = 3 than at 2: = 0, even 
though their mean luminosity at 2 = 3 is ~30 to 40 times larger than that at 
2=0 (Elvis 1994). The short lifetime scenario allows for both density and 
luminosity evolution of the quasar population. However, the ROSAT results 
(Boyle, et al. 1993) appear to be consistent with the simple assumption 
of pure luminosity evolution, and we will use that analysis accordingly in 
updating our calculation of the neutrino background from quasars. 

3. Cosmic Ray Interactions in Quasar Cores 

Let us assume that cosmic rays are accelerated in the centers of active galax- 
ies (see section 4) and examine their expected interactions and energy losses, 
including those which lead to 7-ray and neutrino production. First, we con- 
sider the photon environment. 

A characteristic quasar spectrum has comparable amounts of power in 
each logarithmic energy interval, but the spectrum is not completely flat. 
Typically, the most prominant feature of the spectrum is the "blue bump" , 
which has a pronounced peak in the hard UV and is thought to be thermal 
emission from the black hole accretion disk. This feature is thought to arise 
from the release of energy through viscous processes in the accretion disk 
{e.g., Laor 1990). 

There is also usually a thermal feature in the IR region, thought to be 
due to dust intercepting and reradiating the UV emission. This hypothesis is 
supported by the fact that the observable "blue bump" is not as pronounced 
in AGN which have prominant thermal IR features. Dust can only exist 
outside of the nucleus, as it is destroyed by sputtering, so the IR spectrum 
must arise at a large distance from the "central engine" (Barvainis 1987). 
The lack of variability in the IR spectrum also implies that the emission is 
from a large region (Edelson and Malkan 1987). 

Between the IR and UV emission features there may also be a further 
continuum component (Loska and Czerny 1990; Laor 1990), but this com- 
ponent appears to be not associated with the very central core (Done et al. 
1990). 
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The most rapidly variable part of the spectrum is the X-ray region. Its 
emission is similar in luminosity to that of UV bump but its rapid time 
variability indicates that the X-rays are produced in a smaller region than 
the UV flux (Lawrence, et al. 1987; McHardy and Czerny 1987; McHardy 
1988). Thus, in the very central regions, only the UV and X-ray photon fields 
are of importance. 

The lack of strong X-ray absorption features in quasar spectra (Mushotzky 
1982; Turner and Pounds 1989) implies that the secondary X-rays are pro- 
duced in regions of low column density. This puts strict limits on the amount 
of target gas for pp interactions, and one flnds that the very large photon 
density in the quasar core makes prf interactions the dominant energy loss 
mechanism. This is a threshhold interaction, and the dominant v produc- 
tion channel at energies just above the threshold is 7p — ?■ A — ?■ rayr"*" (Stecker 



The optical depth to photomeson production for protons of energy E,p 
can be approximated by T{Ep) ~ tn{t)aoR where Gq at the A resonance 
peak is ~ 5 X 10~^® cm^ and ra(e) is the differential number density of 
photons at the resonance energy in cm~"^ eV~^. An estimate for ra(e) from 
the accretion disk can be found by assuming that all sources run at some 
constant efficiency (see Section 2). Fitting accretion disk models to the UV 
spectrum (Laor 1990) gives a typical luminosity of (0.03 — 0.1)Ledd ergs s~^. 
We have assumed a "blue bump" UV luminosity of ~ O.OSLgdd is emitted 
isotropically. Given the weak dependence of the characteristic accretion disk 
temperature on black hole mass and the very similar UV bump shapes seen 
in the observational data, we can assume a luminosity independent 'generic' 
UV quasar spectrum. The quasar X-ray spectrum is also fairly universal, 
typically being an e"^''' power law above 2 keV (Mushotzky 1982; Turner 
and Pounds 1989), with a cutoff in the MeV energy range (Rothschild, et 
al. 1983). The total X-ray luminosity is roughly the same as that in the UV 
bump i.e. ~ Luv ~ O.OSLgdd- Normalizing so that L ~ AttR^c J tn{t)dt = 
O.lLeddj we can calculate the mean UV and X-ray photon density within 
the UV emission region of radius R. We thus obtain 



where L45 is the total UV luminosity in units of 10"*^ ergs s~^. The 1/L 
scaling of ra(e) follows from the linear scaling of both R and L with M 
(i.e. with Ledd and R^ above.) 

The optical depth from the "blue bump" UV photons is ^ 1. The optical 
depth is independent of the luminosity since T{Ep) oc n{t)R and ra(e) cx 1/L 
while R oc L. The high optical depth of the UV photon held implies that the 



1968). 




0.165e e < 1 eV 

0.165e-°-9 1< e < 40 eV 

5.35 X lO-^e^e-^/i^ 40 < e < 192 eV 

4.15 X lO-^e-i-^ 192 < e < 10^ eV 



(1) 
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secondary neutrons also will not generally escape the luminous core region, 
so that comparable amounts of power are generated via nj collisions. With 
roughly half of the energy loss going into 7r°'s and the other half into tt^'s, 
the luminosity for z/^ + is ~ OAL^, and for + Fe is ~ 0.2Lx. As the 
quasar photon spectrum drops exponentially between 40 and 200 eV, it is 
overwhelmingly the UV bump photons which interact with the high-energy 
nucleons. We note that for this typical quasar spectrum, the X-rays are not 
optically thick. 

4. Proton Acceleration in Quasar Accretion Disk Shocks 

It has been suggested by Kazanas and Ellison (1986) that protons can be 
accelerated to very high energies by diffusion across the strong shock gen- 
erated at the inner edge of the black hole accretion disk. First-order Fermi 
acceleration of protons in strong (non-relativistic) shocks produces a power- 
law proton energy spectrum cx E~'^ up to a maximum energy i^max- This 
maximum proton energy is determined by equating the proton acceleration 
time, ts,cc{Ep), with the pj lifetime, tp^{Ep). The proton acceleration time is 
given by iacc(-E'p) ~ 2.2 X 10~'^{Rshock/ Rs){Ep/ nip) (Kazanas and Elli- 
son (1986). The pj lifetime is given by tp^{Ep) ~ {N^ap^cK)~^ , where k 
is the mean elasticity and is the number density of photons above the 
photoproduction threshhold. The shock radius, i?shockj is fixed at ~ lOi?^ by 
our assumption that the X-ray luminosity is O.OSLgdd- The magnetic field 
B is taken to be 10"^ L45 G, which assumes approximate equipartition 
of the magnetic energy density with the quasar accretion disk "UV bump" 
photon density. (Note that this is a change from our earlier work (Stecker, 
et al. 1991) where we just assumed a luminosity independent value of B 
= 10"^ G.) Setting the loss time equal to the acceleration time, and noting 
that (X 1/L (see section 3), we obtain -Emax given by equation (2). For 
L45 > 0.1, -Emax cx L^/^, a result similar to that obtained by Szabo and 
Protheroe (1994) using this same assumption about the B field. 

r 0.25L<lf^ L45 < 1 
^max = lO^GeV • <^ O.25LO55, 1< L45 < (2) 
[ 25, L45 > 10^ 

For all luminosities, the maximum proton energy is limited by pj inter- 
actions rather than by the scale size of the shock region. 

5. The Production of Ultrahigh Energy Neutrinos 

Given a power law (roughly cx E~'^) proton spectrum, the pj pion producing 
interactions occur mostly close to the threshhold energy. The dominant v 
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production channel at energies just above the threshold is 7p — ?■ A — ?■ nir'^ 
(Stecker 1968). The cross-section for this reaction peaks at photon energies 
of e' = O.SbnipC^ in the proton rest frame. In the observers frame eEp = 0.35 
eV-EeV with Ep in EeV and e in eV. For UV bump photons, with a mean 
energy of ~30 eV, this translates to a characteristic proton energy of Ec ~ 
10^ GeV. 

It is possible to estimate the isotropic X-ray flux from the core of the clos- 
est and brightest quasar, 3C273. Its X-ray spectrum is observed to exhibit 
an iron line at 6.4 keV of equivalent width ~ 50 eV (Turner, et al. 1990). 
The currently accepted explanation for this line in quasars is fluorescence 
from the X-ray illuminated accretion disk , so that much of the X-radiation 
is probably unbeamed. The isotropic X-ray emission in the 2-10 keV band 
is then ~ 10"*^ ergs s~^, giving a total X-ray flux of ~ 10'*'' ergs s~*. 

It now appears that, in the case of the Seyfert galaxy NGC 4151, at 
least some of the X-ray flux and perhaps most or all of the X-ray flux may 
be of thermal origin with Comptonization playing a major role (Zdziarski, 
Lightman and Maciolek-Niedzwiecki 1993; Titarchuk 1994). Such a contri- 
bution to a quasar X-ray spectrum may correspondingly lower our estimate 
of its neutrino flux, since the neutrino flux in our model is proportional only 
to the non-thermal AGN component of the X-ray emission. (See, however, 
Zdziarski et al., 1994; 1995.) To take account of this eventuality, we will 
assume that, on average, only 30% of quasar X-ray emission is non-thermal 
in normalizing our calculations. Using this estimate, we calculated the neu- 
trino spectrum of 3C273. The result shown in Figure 1 was calculated using 
Monte Carlo techniques. 

The mean pion energy is < i?^ >~ 0.2Ep, while the mean neutrino energy 
is ~ Et^/4 ~ O.ObEp (Stecker 1979). Thus the neutrino spectrum from proton 
interactions with the UV bump extends from Ec/20 to £'inax/20 i.e. from 
~ 10^ to ~ 10® GeV, with a spectral index of ~ 2, mirroring that of the 
proton spectrum. Above ~ 10^ GeV, the cutoff in the proton spectrum leads 
to a cutoff in the neutrino spectrum. Below Ec/20 relativistic kinematics give 
rise to a flat v energy distribution (Stecker 1979). 

It is probable that there is a range of X-ray compactness in quasars, and 
hence a range of optical depths to the pj process. For more compact quasars, 
the maximum energy of the protons would be lower, with the pj interactions 
with X-rays determining the cutoff (Sikora, et al. 1987). Interactions with 
X-rays will also result in the more copious production of lower energy neu- 
trinos with a softer low energy spectrum. This would occur either in models 
with higher compactness (Sikora and Begelman 1992) or indefinitely large 
nucleon containment times (Szabo and Protheroe 1994). The true situation 
is probably somewhat intermediate between the various models. 

A soft spectrum of low energy neutrinos can also arise from the decay of 
pions produced in pp interactions. Using the X-ray data, we have argued that 
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Fig. 1. The predicted i^^ + i^^ flux from 3C273, with Lx = 10^^ ergs s ^ (Piccinotti, et al. 
1982) and redshift z = 0.158. (A Hubble constant of 50 km-s"^-Mpc"^ has been assumed 
in our calculations.) Note that the i^e + i^e flux is half that of the i^^^ + flux. 

such interactions should not account for the bulk of the neutrino production 
(see section 3). In our initial work (Stecker, et al. 1991), we estimated that, at 
most, 10% of AGN neutrinos would be from pp interactions. This conclusion 
is consistent with the suggested accretion disk production estimate of Nellen, 
et al. (1993). 

6. The Extragalactic Neutrino Background Spectrum 

The diffuse v flux can be found by integrating over the luminosity function 
(LF) for quasars. The most recent determinations of the X-ray LF, as well as 
high redshift luminosity evolution, come from the ROSAT satellite (Boyle, 
et al. 1993). We will adopt model I of Boyle, et al. (1993). For this model, 
the LF at 2: = is given by 



where the luminosity function p is per luminosity interval in units of (10'*'* 
erg s~*)~*, L44 is the 0.3 to 3.5 keV X-ray luminosity in units of 10'*'* erg 




(3) 
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s~^, Linin = lO"^^ erg , Li = W^^-^ ergs and Lmax = 10'*'' ergs s~*. 
These luminosities refer to the X-ray flux in the ROSAT range and must 
be multiplied by a factor 17 to transform them into a total X-ray flux. The 
source luminosity as a function of redshift, z, is taken to be (Boyle, et al. 
1993) 

^^^^ ^^(0)(l + ^)"^ Z<2 



We calculate the diffuse neutrino background assuming an epoch of forma- 
tion of quasar activity corresponding to z^^x = 4 and note that our results 
are not very sensitive to z^^x in the range between 2 and 4. 

Figure 2 shows our new background result, obtained assuming that ~30 
% of the quasar X-ray flux is from the jrf interactions and subsequent elec- 
tromagnetic cascading (see section 5). Our results indicate that the diffuse 
background neutrino flux from quasars could dominate over the atmospheric 
background in the range 10^ to 10*° GeV. 



7. Possible Origin of the Broad Line Emission Region 

The broad line region (BLR), a characteristic feature of AGN, is thought 
to be made up of a large number of dense, illuminated clouds (Davidson 
and Netzer 1979). The large line widths of ~ 5000 km s~* can be explained 
by doppler shifts due to orbital motions. If the clouds have predominantly 
circular orbits, then the velocity can be calculated as u = \/GM/r, or /3 = 
\/Rs/1r. For these high velocities, the clouds have to be within ~ 1700i?s, 
or 5 X 10*^ cm for a 10® solar mass AGN black hole. 

An outstanding problem in this scenario is the formation and stability of 
these clouds. While the two phase instability, in which dense condensations 
can exist in pressure equilibrium with a hot, rarifled gas is seen as the best 
and currently accepted model (Krolik, McKee and Tarter 1981), there are 
many problems with this explanation (Mathews 1986; Elizur and Ferland 
1986). 

However, we have previously showed that the high energy neutrino lumi- 
nosities predicted from quasars will have a profound effect on stars close to 
the center of the host galaxy, producing stellar winds, swelling the atmo- 
spheres and even causing their total disruption (Stecker, et al. 1985; Gaiss- 
er, et al. 1986; Stecker, et al. 1991). The column density for absorption of 
a ~ 10^ GeV neutrino is A = vanjoy^ ~ 2 X 10^ g cm~^, while total 
column density for a solar mass star < p > R ^ 10*^ g cm~^. Using a 
conservative disruption criterion that the neutrino energy deposited in the 
star be greater than its nuclear energy generation, we flnd a sphere of stellar 
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disruption, i?ssDj within a radius of 

i?ssD~30L^/'(^)"''' light days (5) 

where M* and Mq are the masses of the irradiated star and the Sun respec- 
tively. This radius is the same as that inferred above for the BLR, and 
matches well with the current best observational determinations (Netzer 
1989; Clavel, et al. 1991) and previous calculations have shown that out- 
flowing material from stellar disruption provides, at least qualitatively, an 
environment similar to that in which the broad lines must arise Kazanas 
1989). Thus, high energy neutrinos from an AGN core could play a key role 
in producing the BLR. 

8. Neutrinos from the Cores and Jets of Blazars 

The evidence for relativistic bulk motions and beaming in blazars ("optically 
violent variable" quasars and BL Lac objects) implies that the intrinsic 
luminosities and sizes of the emission regions are more uncertain than in 
radio quiet objects. However, if the above discussion of the origin of the 
BLR is correct, then the cores of these objects should similarly produce 
neutrino fluxes. These objects would be expected to contain the same sort 
of nucleus as the radio quiet objects, but with the addition of jet emission 
from further out. 

With the observation of 7-ray emission above 100 MeV from at least 40 
blazars out to redshifts greater than 2 by the EGRET group (Fichtel 1994), 
it has now become possible to estimate the corresponding high energy neu- 
trino emission from the jets of these sources. Let us assume that the observed 
7-rays are produced in the relativistic jets of these objects and are beamed 
because (a) as discussed before, we do not expect 7-rays above a few MeV 
to escape from the core of a typical quasar, and (b) EGRET observed time 
variability on the order of a day from 3C279, implying relativistic beaming. 
Let us further assume the the observed 7-rays are produced by high energy 
protons, leading to pion production and the subsequent decay of neutral 
pions into 7-rays and accompanying charged pions into electrons, positrons 
and neutrinos. We further note that the observed differential spectrum of 
3C279 is consistent with a roughly E~'^ power law up to energies in excess 
of 5 GeV (Hartman, et al. 1992), as would be expected for shock accelerated 
proton primaries leading to pion decay 7-rays. Also, the detection of 7-rays 
above 0.5 TeV from Mrk 421 by the Whipple Observatory group (Punch, 
et al. 1992), in combination with the EGRET detection, implies a rough- 
ly E~'^ spectrum for this source up to TeV energies and possibly beyond. 
Such high energies are hard to account for from electron acceleration alone, 
since radiative energy losses of electrons at such energies are quite severe. 
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Mannhiem, et al. (1992) have proposed a model for proton acceleration in 
blazar jets. 

If we make the above assumptions, we expect that the ratio of 7-rays 
to + T'i_t should be 1:1. Assuming an E~'^ differential spectrum, the ratio 
above a given energy should be 2:1, since the mean energy of a 7-ray from 
pion decay is roughly twice the energy of a neutrino from the decay of a 
pion of equal energy. We can then predict TeV neutrino fluxes from the 
EGRET sources. For Mrk 421, with a reported 7-ray flux above 0.5 TeV of 
~ 1.2 X 10~^^ cm~^ s~^, we obtain an immediate predicted neutrino flux 
above 1 TeV of ~ 6x 10~^^ in these same units. This source is at a low enough 
redshift (0.031) so that no significant 7-ray absorption at TeV energies is 
expected, however, the other EGRET sources are at higher redshifts where 
pair production interactions with intergalactic infrared photons can cut off 
the observed 7-ray spectrum (Stecker, De Jager and Salamon 1992). Indeed, 
although Mrk 421 is the weakest of the EGRET sources yet reported, none of 
the other sources searched for were detected by the Whipple team above 0.5 
TeV (Kerrick, et al. 1993). However, correspondingly high energy neutrinos 
can reach us unattenuated, so that we may use an extrapolation of the 
EGRET spectrum to predict TeV neutrino fiuxes. 

For the jet in 3C279, in its more quiescent state, corresponding to the 
second EGRET observing period, extrapolation of the E~'^ 7-1'ay spectrum 
leads to a predicted jet muon neutrino fiux of ~ 4 X 10~^^ above 1 TeV. 
For the fiaring phase the corresponding neutrino fiux would be ~ 2 X 10~^°, 
twice the minimum detectable fiux for DUMAND II. Our predicted neutrino 
fiux above 1 TeV for the core of this source would be ~ 10"^"^. 

For the 3C273 jet, extrapolating the EGRET spectrum (C. von Montigny, 
et al. 1993) leads to a predicted z/^ fiux above 1 TeV of ~ 10"^"^. The 
corresponding prediction for the core fiux is ~ 4 X 10"^"*. 

It is interesting to note that for superluminal jet sources, with a typical 
bulk Lorentz factor of ~ 10, the beamed fiux would appear to be ~ 10"^ times 
brighter than that of an isotropic source with the same intrinsic production 
rate. 

Using the 1:1 correspondence of 7-rays to z/^'s (see above), combined with 
an extrapolated estimate of the blazar 7-ray background which we published 
recently (Salamon and Stecker 1994), we immediately obtain the estimate of 
the neutrino background fiux from blazar jets shown in Figure 2, assuming 
the spectrum continues as a power law well above TeV energies. Figure 2 
also shows the AGN neutrino background fiuxes calculated by us here and 
by other workers. 



agnnu.tex - Date: January 17, 1995 Time: 11:07 



12 



F.W. STECKER AND M.H. SALAMON 




3456789 10 
logE,(GeV) 



Fig. 2. Our revised integrated high energy i^^+F^^ neutrino background from quasars (Q) 
and blazars (B) (thick soUd Unes), as well as the fluxes calculated by Sikora and Begelman 
(f992) (dotted line), Biermann (f992) (dashed line), and the geometric mean for the h = I 
model of Szabo and Protheroe (f994) (dash-dot line). Also shown is the horizontal i^^+F^^ 
flux from high energy cosmic rays interacting with the Earth's atmosphere (ATM, solid 
line) (Stecker f979). 

9. Discussion 

As can be seen from Figure 2, the high energy neutrino background fluxes 
from AGN calculated by various workers are comparable. This is not too 
surprising, since they are mainly based on the assumptions of shock accel- 
eration of protons to high energies followed by photomeson production and 
subsequent charged pion decay. All results are normalized to AGN X-ray or 
X-ray background energy fluxes in some way. 

The results shown in Figure 2 lead to predictions of observable high ener- 
gy neutrino background fluxes. Indeed, presently existing data are beginning 
to test the different predictions of the models (Zas, Halzen and Vazquez 
1993; Rhode 1994; Aglietta, et al. 1994; Gaisser, Halzen and Stanev 1995). 
In particular, the analyses of upper limits on muon fluxes from the Frejus 
experiment (Rhode 1994; Gaisser, et al. 1995) are consistent with our quasar 



agnnu.tex - Date: January 17, 1995 Time: 11:07 



HIGH ENERGY NEUTRINOS FROM QUASARS 



13 



background flux prediction, but appear to rule out part of the range of mod- 
el parameters considered by Szabo and Protheroe (1994). Our prediction is 
also consistent with bounds on giant horizontal air showers established with 
the AKENO array (Zas, Halzen and Vazquez 1993) and the EAS-TOP col- 
laboration (Aglietta, et al. 1994). However, the predicted flux is not very far 
below the presently established bounds. 

The event rates for various v detectors can be calculated from our back- 
ground spectrum given in Figure 2. The Fly's Eye detector, with its energy 
threshhold of 10® GeV would expect to see ^ 1 event per year, consis- 
tent with their present non-detection. Their next generation instrument, 
the HiRes Eye would see < 1 downward events per year. Above the HiRes 
energy threshold of ~ 10® GeV, upward neutrinos are severely attenuated by 
absorption in the Earth. However, the upward event rates at energies greater 
than 10^ GeV for DUMAND H, AMANDA and NESTOR are of order 10^ 
yr-^ 

We note that n^f interactions produce a signiflcant flux of z7e's at 6.3 PeV, 
the Glashow resonance energy for the Vee~ — ?■ W~ process. An observation 
of an enhanced event rate at this energy would indicate that AGN are opti- 
cally thick to neutrons. We estimate that DUMAND H would see on the 
order of 10 events per year from W~ resonance production. (The reader is 
referred to the articles in Stenger, et al. (1992) for a discussion of various 
neutrino detector experiments.) 

We conclude that results from high energy neutrino detectors either 
planned or presently under construction should be able to test the viabili- 
ty of the idea that signiflcant proton acceleration in AGN cores ultimately 
produces a substantial fraction of their observed X-ray luminosity. 
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